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ABSTRACT 


The  current  aeter  records  collected  at  three  sites  in 
the  Gulf  of  Mexico  during  the  passaqe  of  Hurricane  Prederic 
are  analyzed  to  deteraine  the  s tora- induced  flow  at  various 
ocean  depths,  deteraine  the  associated  energy  increase  and 
decay,  and  compare  these  observations  to  siailar  results 
froa  a  numerical  model.  The  records  at  the  two  deeper  sites 
are  rather  unique  because  they  are  within  100  km  of  the  hur¬ 
ricane  track.  Pre-storm  conditions  are  controlled  by  topog¬ 
raphy,  and  as  the  stora  passes  there  is  an  abrupt  change  in 
the  direction  of  flow  and  initiation  of  a  strong  inertial 
response  at  all  levels  of  the  two  deeper  sites.  After  this 
initial  surge,  the  residual  flow  tends  toward  the  pre-storm 
direction.  The  horizontal  kinetic  energy  associated  with 
inertial  motion  is  calculated.  The  energy  increase  and 
decay  is  shown  to  vary  with  depth. 

An  embedded  mixed-layer  ocean  circulation  model  (Adamec 
et  al#  1981)  is  forced  with  an  idealized  storm  translating 
at  the  same  speed  (7.5m  s~l)  as  Prederic.  The  abrupt 


response  and  strong  inertial  component  predicted  by 
model  is  qualitatively  siailar  to  the  observations. 
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I.  IHTBOp  OCTION 


Several  studies  of  the  response  of  the  ocean  to  hurri¬ 
cane  passage  have  been  made  by  such  authors  as  Fisher 
(1958),  Leipper  (1967),  Bright  (1969),  Pudov  et  §1  (1979) 
and  Fedorov  gt  ai  (1979).  A  review  of  the  hydrographic, 
^.e.  sea-surface  temperature  (SST)  ,  salinity  tenperature 
depth  (STD),  and  expendable  b  athytheraograph  (XBT)  ,  sur¬ 
veys  of  several  of  these  authors  is  given  in  Table  I  froa 
Price  (1981).  The  aajority  of  these  observations  concen¬ 
trate  on  the  asymmetrical  SST  response,  and  the  expected 
ocean  mixed-layer  response  to  hurricane  passage  of  vertical 
mixing  as  the  hurricane  approaches,  followed  by  upwelling  in 
the  wake  of  the  storm  (see  e.^.  Leipper,  1967;  Friase, 
1977).  Additionally,  oscillations  on  the  order  of  the  local 
inertial  period  have  been  noted  in  the  temperature  and 
current  fields  in  the  wake  of  open-ocean  storms.  The  iner¬ 
tial  response  becomes  complicated  in  the  shoal  waters  of 
continental  shelves  and  the  ability  to  detect  inertial 
moticn  is  severely  restricted  (Hayer  et  al,  1981).  The  hur¬ 
ricane-induced  upwelling  decreases  the  mixed  layer  depth. 
The  depth  reaches  a  local  minimum  about  one  half  inertial 
period  after  eye  passage  and  continues  to  oscillate  at  about 
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the  local  inertial  period  for  some  period  after  the  storm 
passage.  Geisler  (1970)  presented  the  relationship  between 
hurricane  translation  speed  and  the  oscillatory  response  of 
the  theraocline.  If  the  hurricane  translation  speed  is 
greater  than  the  internal  phase  speed,  typically  about  2  o 
s”1,  and  if  its  horizontal  scale  is  comparable  to  the  inter¬ 
nal  Kossby  radius,  currents  throughout  the  affected  area  are 
controlled  by  a  balance  between  the  centrifugal  and  Coriolis 
accelerations. 


TABLE  I 

Hydrographic  Studies  of. the  Sea  Syrf ace^ Temperature  Response 

to  Hurricanes  (Price,  1991): 


Humean* 


Study:  Humean* 

Method:  Region 

Average 

U. 

(■  a'1) 

central 

pnaaura 

'  (mb) 

ASST.. 

CO 

Position  of*  A  SST^ 

Leippcr  ( 1967):  HiMa(1964) 

extensive  post- hurricane 
hydrographic  survey:  Gulf  of 
Mexico 

} 

930 

-6 

pattern  ia  generally 
unclear,  may  be  50  has 
to  left  of  track  (Fig.  1) 

Fedorov  ,i  a/.  (1979): 

Ella  (I960 

Extensive  pro-  and  post-hurricane 
XBT  survey:  mid-Atlantic 

6 

9(0 

-2 

30  km  to  right  (Fig.  3. 
same  as  this  Fig.  la) 

Pudov  ft  al.  (1979):  Ttu 
(I97J) 

extensive  post- hurricane  STD  survey: 
mid-pacific 

6 

940 

-4 

75  km  to  right  (Fig.  1. 
same  as  this  Fig.  2a) 

Wright  <  1969):  Shirlty  ( 1965) 

1  pre-,  1  pott-hurricanc  XBT  section: 
vicinity  of  the  Kuroshio 

15 

935 

-3 

20  km  to  right  (Fig.  4) 

Jordan  (1964): 

extensive  pre-  and  post-hurricane 

16 

920 

-2 

130  km  to  right 

Wanda  (19561 

Clara  (1955) 

SST  reports  from  ships  of 
opportunity:  mid- Pacific 

IS 

915 

50  km  to  right 
(Figs.  2  and  3) 

^Estimates  made  by  Price  (19  81)  from  their  figures  noted 


the 


The  response  to  Hurricane  Eloise  as  her  eye  passed  over 
National  Oceanic  and  Atmospheric  Administration’s  Data 
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Buoy  Office  (NDBO)  EB-10  buoy  has  been  studied  extensively 


(see  £.q.  Martin,  1982;  Price,  1981;  Black  and  Bithee, 
1976).  EB-10  was  located  in  the  central  Gulf  of  Mexico  and 
collected  the  first  open-ocean  data  under  a  hurricane  as 
reported  by  Bithee  and  Johnson  (1976)  .  Mayer  et  a*,  (1981) 
reported  a  study  of  the  passage  of  Hurricane  Belle  over  the 
continental  shelf  of  the  New  York  Bight  on  10  August  1976. 
These  continental  shelf  responses  differ  from  the  open-ocean 
respcnses  due  to  the  topographic  influences  of  the  shelf  and 
large  gradients  in  the  physical  properties  between  the  shelf 
water  and  the  deep  ocean  water. 

Hurricane  Frederic  passed  through  the  Gulf  of  Mexico 
between  11  and  13  September,  1979.  The  satellite  data 
depicted  in  Figs.  1,  2,  and  3  together  with  whatever  XBT 
data  collected  during  hurricane  passage  are  usually  the  only 
type  of  data  available  to  study  these  geophysically, 
socially  and  economically  important  events.  However,  during 
this  period  the  Naval  Oceanographic  Office  (NAVOCEANO)  ,  NSTL 
Station,  Bay  St.  Louis,  Mississippi  operated  three  data 
buoys  which  collected  a  unique  set  of  ocean  temperature  and 
current  data.  This  thesis  studies  the  ocean  current 


response  to  the  passage  of  Hurricane  Frederic  as  it  passed 


near  the  three  data  buoys.  Additionally  a  comparison  of  the 
real  data  to  the  results  provided  by  the  three-dimensional 
ocean  model  of  Adaaec  si  il  (1981)  is  made.  The  hurricane 
forcing  of  the  model  is  idealized  rather  than  being  formu¬ 
lated  to  represent  the  actual  hurricane,  but  a  translation 
speed  of  7.5  a  s-*,  the  same  as  for  Frederic,  is  used. 

The  buoy  current  data  provided  by  NAVOCEANO  indicate 
that  the  floe  associated  with  the  passage  of  Hurricane 
Frederic  had  a  large  inertial  component  after  the  storm 
passed.  The  rate  at  which  the  inertial  flow  damps,  and  the 
rate  at  which  it  propagates  with  depth,  are  determined  from 
the  buoy  data  and  compared  to  the  model. 
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II.  DATA  PROC3DOR  SS  4£D  aEgOJjIS 


A.  EACKGROUND 

Ten  current  meter  array  records  collected  at  three  sites 
on  the  continental  shelf  of  the  southeastern  United  States 
(Pig.  4)  during  the  passage  of  Hurricane  Prederic  were  pro¬ 
vided  by  NAVOCEAUO.  As  depicted  in  Fig.  5,  the  arrays  con¬ 
sisted  of  Aanderaa  spar  buoys  anchored  to  the  bottom.  The 
anchor  cable  held  instruments  at  three  depths  at  each  of  the 
stations  one  and  two  and  at  four  depths  at  station  three. 
The  shallowest  instrument  depth  was  19  m  at  station  two, 
while  the  deepest  was  457  a  at  station  three.  The  bathyme¬ 
try  (Fig.  5)  shows  all  three  stations  are  or.  the  continental 
shelf.  Data  were  provided  for  a  continuous  period  from  a 
few  days  before  the  storm  passage,  about  2100  Greenwich  Mean 
Time  (GMT)  12  September  1979,  until  a  few  days  after  the 
storm  passage. 

Water  depth  at  station  one  is  only  about  100  meters. 
Price  (1981)  reported  difficulties  with  meters  in  a  similar 
coastal  environment  because  of  the  presence  of  very  strong 
horizontal  gradients  and  strong  topographic  effects  on  the 
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Figure  4.  Best  position  track  of  Hurricane  Frederic  and 
positrons  of  current  meter  arrays  (provided  by 
t.K.  shay,  NA70CE&N3)  . 

currents,  (layer  et  al  reported  intense,  first-mode,  near- 
inertial  period  oscillations  in  water  of  about  70  meter 
depth.  Similar  observations  at  sites  in  shallower  water  of 
about  50  meter  depth  resulted  in  only  weak,  heavily  damped, 
second-mode  oscillations. 

The  current  meter  records  used  in  this  thesis  also  show 
a  variation  in  response  among  the  stations  with  different 
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bottom  depths.  At  station  three  where  the  water  depth  is 
about  465  meters,  an  increase  in  current  aagnitude  starting 
about  1920  GMT  12  September  1979  is  very  apparent  in  the  raw 
data  record  (Fig.  6a).  This  record  also  shows  an 
oscillation  which  cycles  through  one  wavelength  about  once 
every  24  hours  (1440  minutes).  This  oscillation  decays 
smoothly  over  about  a  week,  and  at  first  glance  it  appears 
to  have  a  period  of  about  22  hours.  A  more  detailed  analy¬ 
sis  of  the  record  shows  the  zero  down-crossings  to  be  very 
nearly  1450  minutes  apart  during  the  largest  oscillation  in 
the  record.  The  average  period  for  seven  cycles  during  this 
interval  is  1430  minutes  (23.3  hours)  between  down-cross¬ 
ings.  The  local  inertial  period  at  station  three  is  1473 
minutes.  Shortly  after  the  surface  response,  the  rapid 
increase  in  current  magnitude  is  also  obvious  at  the  depth 
of  457  meters  at  station  three  (see  Fig.  6b)  .  The  average 
period  of  the  oscillation  at  this  greater  depth,  1370  min¬ 
utes  or  22.8  hours,  is  less  than  the  period  near  the  sur¬ 
face,  but  the  damping  of  the  oscillation  is  not  obvious  at 
457  meters.  At  station  one  whare  the  bottom  depth  is  only 
about  100  meters,  the  increase  in  current  speed  is  detecta¬ 
ble  in  the  raw  data  record  near  the  surface  (Fig.  7a) ,  but 
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any  period  of  oscillation  on  the  order  of  24  hours  or  a 
damping  rate  is  difficult  to  detect.  An  increase  ir.  current 
velocity  on  12  September  1979  near  the  bottom  at  station  one 
(Pig.  7b)  is  not  apparent  and  it  can  be  seen  that  current 
magnitudes  are  larger  on  9  September  than  they  are  on  12 
September. 

To  determine  if  the  observed  currents  were  extraordi¬ 
nary,  some  idea  of  the  likely  currents  in  the  area  of  obser¬ 
vations  is  helpful.  A  dominant  feature  of  the  3ulf  of 
Mexico  surface  currents  is  the  Loop  Current  (Pig.  8).  This 
is  a  clockwise  current  of  about  50  to  203  cm  s~l  and  90  to 
150  km  width  (Leipper,  1970)  which  enters  the  Gulf  in  the 
west  as  the  Yucatan  Current,  and  exits  through  the  Plorida 
Straits  as  part  of  the  Gulf  Stream  system.  The  Yucatan 
Current  flows  north  from  Honduras  between  the  Yucatan 
Peninsula  and  Cuba  into  the  central  east  Gulf  and  forms  the 
western  section  of  the  Loop.  Further  north  the  current 
flows  east,  southeast,  and  east  again.  Elsewhere  in  the 
northeastern  Gulf,  the  surface  flow  is  generally  cyclonic 
along  the  coast  until  turning  southward  near  Louisiana. 
Eddies  are  known  to  have  detached  from  the  Loop  Current 
(Leipper,  1970)  and  could  move  into  the  area  of  the 
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instrument  arrays.  Leipper  (1970)  shows  a  systematic  devel¬ 
opment  and  breakdown  of  the  Loop  Current  which  places  the 
northern  edge  of  the  current  further  north  in  the  spring 
than  in  other  seasons  (Fig.  9).  Holinari  (1978)  suggests 
the  northward  intrusion  of  the  Loop  Current  is  not  seasonal, 
but  that  climatological  results  are  biased  by  temporal  sam¬ 
pling  techniques.  Hurlburt  and  Thompson  (1980)  presented  a 
numerical  study  of  Loop  Current  intrusions  and  eddy  shed¬ 
ding.  They  obtained  theoretical  expressions  for  the  eddy 
diameter  and  penetration  distance  of  the  Loop  Current  into 
the  Gulf  and  concluded  that  if  only  vorticity  dynamics  are 
considered,  the  interaction  between  the  Florida  Shelf  topog¬ 
raphy  and  the  pressure  field  results  in  a  balance  which 
stops  the  northward  penetration  of  the  Loop  current.  All  of 
the  current  array  stations  are  further  north  than  29  deg 
North,  and  although  the  waters  of  the  Loop  rarely  extend 
this  far  north.  Huh  at  al  (1981)  examined  an  intrusion  of 
Loop  waters  as  far  north  as  the  data  buoys  which  collected 
the  raw  data  used  in  this  thesis.  It  is  not  possible  to 
tell  a  priori  whether  or  not  the  Loop  current  is  in  the  area 
of  the  current  meter  arrays  during  the  period  of 
observations,  but  it  seems  unlikely. 
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Figure  8. 
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Another  factor  which  might  be  expected  is  the  influence 


of  topography.  The  principle  cf  conservation  of  potential 
vorticity  causes  the  current  to  follow  the  bathymetric  con¬ 
tours.  The  proximity  cf  the  buoys  to  the  coast  suggests 
that  near-shore  property  gradients  between  the  shelf  water 
and  the  deep  water  may  have  important  dynamic  and  thermody¬ 


namic  consequences.  The  Desoto  Canyon  is  also  in  the  area 


cf  observations 


and  may  tend  to 


enhance  cross-shelf 


Figure  9.  The  spring  intrusion  of  1966  as  indicated  by 

overlays  of  the  150-aeter  contour  lines  from  the 
22°C  topographies  of  all  spring  cruises  of  1966. 
Hay  1964  and  November  and  December  1965  ape 
included  for  supplemental  information  (Leipper, 
1970). 


B.  HAW  DATA 


The  NAVOCEANO  data  consisted  of  north-south  (v)  and 


east-vest  (a)  velocity  components  and  temperature  (X)  at  10 


minute  intervals.  Temperatures  above  21.5  °C  were  not 


recorded  due  to  thermistor  limitations.  Temperatures  above 


21.5  °c  occurred  only  near  the  surface  and  complete  tempera¬ 
ture  records  were  available  only  from  the  deeper 
instruments.  Also,  NAT  OCEANS  reported  apparent  internal 
clock  synchronization  problems  in  the  current  meter  record 
at  64  meters  depth,  station  one. 

The  hurricane  passed  the  data  buoys  with  a  translation 
speed  of  about  7.5  m  s_l,  and  as  depicted  in  Pig.  4,  the  eye 
passed  about  100  km  to  the  west  of  station  three  at  2200  GMT 
12  September  1979  (Julian  day  2  55).  Maximum  winds  of  115 
kts  occurred  at  a  radius  of  about  30  km  from  the  storm  cen¬ 
ter.  The  inertial  periods  for  the  data  buoy  locations  range 
from  23  hours,  57  minutes  at  station  one  to  24  hours,  33 
minutes  at  station  three.  That  is,  the  inertial  period  at 
all  of  the  stations  is  approximately  one  day.  Since  the 
diurnal  tidal  period  is  24.8  hours  and  the  semi-diurnal 
tidal  period  is  12.4  hours,  difficulties  could  arise  in  sep¬ 
arating  any  inertial  motion  from  tidal  motion  before  the 
storm  arrival.  A  spectral  analysis  would  most  likely  not 
have  sufficient  resolution  to  separate  inertial  and  diurnal 
motion.  The  first  harmonic  of  the  inertial  frequency  would 
also  be  indistinguishable  from  the  semi-diurnal  tidal  fre¬ 
quency.  However,  the  large  increase  in  currents,  and  thus 
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kinetic  energy,  observable  in  the  raw  data  records  would  be 
easily  detected  in  an  energy  spectrum. 

The  initial  data  records  provided  by  NA70CEAHQ  span 
non-coincident  periods.  However,  all  records  contained  data 
froe  1920  GMT  07  September  1979  to  1820  GMT  21  September 
1979.  The  raw  data  records  (Pigs.  6  and  7)  show  an  obvious 
response  to  the  hurricane  forcing  starting  about  12 
September  1979.  This  hurricane  forcing  results  in  increased 
currents,  with  maximum  values  of  130  om  s~l  near  the  sur¬ 
face.  Comparison  of  Pigs.  6a  and  5b  shows  that  the  forcing 
response  was  transmitted  as  deep  as  457  jeters  in  less  than 
one  half  an  inertial  period.  The  response  to  hurricane  pas¬ 
sage  is  more  difficult  to  detect  in  most  of  the  raw  data 
records  of  stations  one  and  two  and  it  is  therefore  not  pos¬ 
sible  to  estimate  the  rate  at  which  the  response  is  verti¬ 
cally  propagated  with  depth  merely  by  looking  at  the  raw 
data  record.  The  ability  to  detect  the  increase  of  energy 
associated  with  storm  passage  at  any  station  and  at  any 
level  is  discussed  in  the  section  titled  ’’Procedures  for 
determining  relative  energies".  Although  this  thesis  does 
not  analyze  the  temperature  records  directly,  the  response 
to  hurricane  passage  is  also  apparent  in  the  temperature 


record  (Fig.  10).  The  temperature  record  shows  an 
oscillation  in  the  thermocline  and  an  apparent  cooling  at 
the  437  a  depth  at  station  three,  followed  by  a  gradual 
waning  to  the  end  of  the  period. 


Figure  10.  Raw  temperature  record  for  7-21  September  1979 
at  a  depth  of  437  m  at  station  3. 

C.  DATA  MANIPULATION 

The  initial  records  were  shortened  to  include  only  the 
period  1920  GMT  7  September  1979  to  1820  GMT  18  September 
1979.  This  period  was  the  longest  possible  period  included 
by  all  records  and  was  a  logica 1  choice  since  the  hurricane 
eye  passed  all  three  buoys  between  2100  GMT  12  September  and 
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0300  GMT  13  September  1979.  By  examining  plots  of  the  data 
(§.g.  Pigs.  6  and  7)  it  was  determined  that  significant 
stora  effects  with  soae  pre-storm  and  soma  post-storm  data 
would  be  included  in  the  shortened  record. 

Progressive  vector  diagrams  (PVD's)  were  plotted  for 
each  record  during  the  above  period.  These  diagrams  indi¬ 
cated  that  the  inertial  motion  was  superposed  on  a  mean 
current  which  was  quite  different  at  each  location.  The 
mean  current  was  not  of  primary  interest  and  attempts  were 
made  to  remove  the  mean  current  from  the  records  so  that  the 
magnitude,  damping  period,  and  associated  energy  of  the 
inertial  motion  could  be  determined.  In  the  first  attempt, 
the  mean  current  over  the  first  three  inertial  periods 
(pre-stora  passage)  was  subtracted  from  all  of  the  u  veloc¬ 
ity  components  and  the  v  velocity  components.  The  three 
inertial  period  average  started  at  1920  GMT  07  September, 
which  is  five  days  before  the  stora  center  passed  the  near¬ 
est  buoy.  Comparison  of  the  initial  PVD's  (Fig.  11)  and 
similar  diagrams  of  the  same  iata  with  the  three  inertial 
period  average  removed  (Pig.  12)  was  made. 

The  initial  PVD's  indicate  a  general  transport  in  direc- 
:ions  which,  in  most  cases,  are  nearly  the  same  before  and 
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are  12.  PVD's  resulting  from  the  removal  of  the  gre- 

storm  averages  over  three  inertial  perioas  from 
the  PVD’s  of  Pig.  11. 


after  the  storm-induced  perturbation.  The  perturbations  are 
generally  reminiscent  of  inertial  motion,  i.§.  oscillations 
with  a  period  equal  to  the  local  inertial  period.  PVD's  of 
purely  inertial  motion  would  transcribe  a  simple  circle  for 
each  inertial  period.  The  circles  would  overlap  if  more 
than  one  inertial  period  were  plotted.  Where  the  first 
three  days  of  a  record  were  similar  to  the  mean  flow,  remov¬ 
ing  the  three  inertial  period  average  from  the  record 
resulted  in  more  nearly  overlapping  circular  PVD's.  Remov¬ 
ing  the  three  inertial  period  average  was  not  effective  in 
isolating  the  inertial  component  of  the  flow  when  the  first 
three  days  of  the  record  were  not  obviously  in  the  general 
transport  direction.  Therefore,  the  mean  current  was  chang¬ 
ing  significantly  during/following  the  storm  at  many  of  the 
current  meter  locations.  &  more  effective  method  of  isolat¬ 
ing  the  inertial  component  from  each  of  the  records  was 
required. 

D.  ISOLATING  THE  INERTIAL  MOTION 

In  a  further  attempt  to  isolate  the  large  inertial  com¬ 
ponent  that  was  obvious  in  the  original  records,  a  running 
average  over  the  inertial  period  was  computed  each  three 
hours  and  then  subtracted  from  the  initial  PVD  records.  The 
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averaging  technique  is  necessarily  applied  at  the  central 
point  of  each  inertial  period,  and  as  a  result  the  half¬ 
periods  at  the  beginning  and  9nd  of  the  records  were  lost. 

However,  the  transient  response  is  not  great  at  these  times 

and  the  information  lost  is  inconsequential.  Next,  the 

average  inertial  period  record  is  linearly  interpolated 

between  3-h  values  to  obtain  the  sane  number  of  data  points 

as  the  initial  PVD  data.  PVD's  of  the  running  inertial 

period  average  data  (e.q.  Fig.  13)  show  that  the  majority  of 

the  inertial  flow  in  the  raw  data  is  removed  by  the  running 

average.  The  running  inertial  period  averaged  PVD's  show 

the  general  transport  direction  as  well  as  some  features  j 

which  should  be  noticeable  in  a  model  of  the  ocean  response. 

&  surge  to  the  northeast,  onshore,  is  noticeable  in  the 
record  of  Fig.  13a,  and  there  is  some  indication  that  the 
post-storm  current  is  greater  than  pre-storm,  i.e.  the  hori¬ 
zontal  displacement  between  inertial  period  marks  is  greater 
after  the  storm  surge  than  before.  The  best  example  of 
these  features  in  any  of  the  records  is  shown  in  Fig.  13b. 

A  surge  to  the  west  followed  by  a  counterclockwise  loop  back 
into  the  pre-storm  flow  direction  is  evident.  Additionally, 
the  post-storm  flow  is  of  much  greater  magnitude  than  tne  I 
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pre-storm  flow.  These  diagrams  were  obviously  a  better  rep¬ 
resentation  of  the  mean  flow  throughout  the  record  than  the 
previously  computed  three  inertial  period  averages.  These 
PVD's  also  show  that  the  mean  current  was  indeed  influenced 
by  the  bottom  topography.  Comparison  of  Fig.  4  and  repre¬ 
sentative  DVD's  at  stations  one,  two  and  three  (Figs.  14  and 
15)  show  that  the  mean  flow  follows  the  bathymetric  con¬ 
tours.  The  current  moves  rapidly  toward  the  northwest  at 
station  one  during  and  for  some  time  following  the  period  of 
hurricane  passage  (Fig.  14a).  A  counterclockwise  loop  which 
would  be  associated  with  inertial  motion  is  not  seen  as  it 
is  in  the  staticn  two  record  (Fig.  14b)  .  Comparison  of 
Figs.  13b  and  15  shows  the  flow  is  along  the  contours  but  in 
opposite  directions  at  the  251  and  the  437  or  457  meter 
depths.  Also,  a  shoreward  surge  similar  to  that  at  the 
near-surface  level  of  station  three  (Fig.  13a)  is  seen  at 
least  one  inertial  period  later  near  the  bottom  (Fig.  15). 
The  tendency  to  return  to  the  pre- storm  flow  direction  is 
much  slower  at  the  near-bottom  records  and  is  not  complete 
by  the  end  of  the  record. 

The  running  inertial-period  averages  ire  then  subtracted 
from  the  initial  PVD  data  point-by-point  and  PVD's  of  the 
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resultant  data  are  plotted  (Pig.  16)  .  Examination  of  these 
resultant  PVD’s  shows  that  station  thrae  exhibits  the  most 
obvious  inertial  response,  possibly  due  to  less  coastal  or 
bottom  influence.  The  resultant  PVD's  for  stations  one  and 
two  (not  shown)  also  show  large  inertial  responses.  The 
station  one  record  still  exhibits  a  strong  flow  which  is 
non-inertial.  The  station  two  record  shows  that  the 
majority  of  the  non-inertial  motion  is  removed. 

It  is  not  possible  from  these  plots  to  determine  pre¬ 
cisely  when  the  response  to  the  hurricane  forcing  first 
begins  at  each  depth.  Since  the  rate  at  which  the  energy  is 
transmitted  vertically  is  of  interest,  the  next  procedure  is 
to  calculate  the  energy  associated  with  each  inertial  period 
at  all  available  levels. 

E.  PROCEDURES  FOR  DETER  HI NING  RELATIVE  ENERGIES 

The  PVD* s  which  were  computed  in  the  above  steps  indi¬ 
cate  the  average  of  the  current  record  is  not  zero.  Also,  a 
breakdown  of  these  records  into  inertial  periods  indicate 
this  nonzero  mean  current  which  does  exist  is  not  constant 
throughout  the  record.  An  interactive  computer  program  was 
developed  to  remove  the  mean  current  from  a  record  consist¬ 
ing  of  inertial  flow  superposed  on  a  mean  current.  The  mean 


velocity  is  calculated  as  the  distance  between  the  starting 
and  ending  points  of  the  record  divided  by  the  duration  of 
the  record.  An  equal  portion  of  this  aean  current  is  sub¬ 
tracted  froa  each  ten  ainute  interval  of  the  record.  The 
aean  current  is  reaoved  from  the  record  one  inertial  period 
at  a  time  rather  than  for  an  entire  record  since  the  mean 
current  is  not  constant  throughout  the  entire  record. 

A  circle  of  radius  8  was  fit  to  each  of  the  inertial 
period  records.  The  circle  coordinates  (xfy)  are  computed 
using  a  mean  velocity.  Given  that 

x  *  R  co  s  (0) 
y  =8  sint£) 

The  aean  velocity  is  coaputed  as  the  time  derivative  of  the 
position  as: 

u  =  dx/dt  =  -R  sin(@)  d0/it 
y  *  dy/dt  =  8  cos(g)  d0/3t 

where 

d@/dt  *  -2  7r/(N  dt) 

and  N  represents  the  number  of  points  used  to  make  the  cir¬ 
cle.  Successive  x  and  y  coordinates  of  the  circle  are  then 
coaputed  as: 

xl  *  xO  ♦  (u  it) 
yl  *  yO  ♦  (v  dt) 
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The  starting  position  is  x  *  0.  and  y  =  - 1  - ,  i.e.  at  270 
degrees  of  a  360  degree  circle  centered  at  (0.,0.),  and  the 
successive  positions  are  created  in  a  clockwise  fashion. 
The  radius  of  the  circle  is  adjusted  using  an  error  mini¬ 
mization  program.  This  program  calculates  the  distances 
between  the  positions  required  to  create  the  circle  and  the 
positions  used  to  compute  the  P7D's,  and  selects  the  circle 
which  gives  the  smallest  sum  of  positional  distances.  One 
fourth  of  the  maximum  dimension  of  the  plotted  data  is  used 
as  the  initial  radius  for  the  minimization  subroutine, 
while  the  radius  is  held  constant,  the  sum  of  the  difference 
between  an  "average  error"  and  an  "individual  position 
error"  is  computed  and  divided  by  the  number  of  differences. 
This  error  difference  is  used  in  the  following  calculations. 
The  distances  between  all  points  of  the  computed  circle  and 
the  data  are  first  compared  to  locate  the  data  point  and  the 
circle  point  which  are  closest.  These  points  are  used  as 
the  starting  positions  for  the  distance  error  computations. 
The  "average  error"  is  the  average  of  the  distances  between 
the  consecutive  points  of  the  data  and  the  computed  circle. 
The  "individual  position  error"  is  the  distance  between  a 
point  on  the  circle  and  a  data  point.  The  error  difference 
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for  cne  radius  is  stored  and,  the  radius  is  slowly  incre¬ 
mented  to  twice  the  original  magnitude,  storing  the  error 
difference  at  each  increment.  The  radius  corresponding  to 
the  minimum  value  of  error  difference  is  used  to  create  the 
circle  which  best  fits  the  plotted  data. 

The  circle  created  using  this  adjusted  radius  is  used  as 
the  best  fit.  It  is  then  possible  to  estimate  the  energy 
per  unit  mass  associated  with  each  inertial  period.  The 
inertial  velocity  (V=fR)  associated  with  the  circle  radius  5 
and  coriolis  parameter  f  is  computed.  One  half  of  the 
square  of  this  velocity  is  assumed  to  be  the  energy  per  unit 
mass.  The  circle  fitting  error  difference  computed  as  above 
is  divided  by  the  radius  of  the  fitted  circle  and  referred 
to  as  the  relative  RHS  error.  The  energy  values  are  plotted 
versus  the  inertial  period  at  each  level  at  each  station 
(See  Figs.  17  and  18).  The  radius  for  each  inertial  period, 
the  relative  RHS  error  of  fitting  the  circle  and  the  mean 
current  removed  from  each  inertial  period  are  listed  in 
Table  II  for  the  records  from  station  one  at  a  depth  of  49  m 
and  from  station  three  at  a  depth  of  21  m.  Also,  the  mini¬ 
mization  program  is  tested  using  exact  circles  as  input 
data.  The  resultant  errors  for  given  radii  are  also  listed 
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in  Table  II.  The  exact  circle  error  shows  the  error  is  near 


zero  when  the  notion  is  circular.  The  error  is  found  to  be 
a  smaller  percentage  of  the  radius,  i. e.  the  relative  RMS 
error  is  snaller,  at  station  three  than  at  station  one. 
This  supports  the  previous  results  of  the  motion  being  more 
nearly  circular,  and  thus  more  inertial,  at  station  three 
than  at  station  one. 

A  good  fit  between  the  positions  of  the  generated  circle 
and  the  positions  from  the  real  data  with  the  inertial  mean 
removed  is  obtained  for  the  records  of  station  three  where 
the  water  is  deepest  and  the  coastal  influence  is  least 
(Fig.  19a)  .  The  agreement  in  positions  between  the  two 
curves  at  station  one  is  generally  less  since  these  records 
are  affected  more  by  the  bottom,  and  by  tidal  motion  which 
is  not  necessarily  removed.  Obviously  non-inertial  motion 
remained  in  the  station  one  record  (Fig.  19b).  Example 
plots  of  the  record  after  all  means  were  removed  are  shown 
in  Fig.  20.  The  figures  show  nearly  overlapping  circular 
motion  except  for  different  radii.  Plots  of  the  energy  val¬ 
uer  versus  inertial  periods  are  given  in  Fig.  17  for  the 
level  nearest  the  surface  at  each  of  the  three  stations. 
Energy  values  of  4000  at  station  3  correspond  to  an 


6 


TABLE  II 


Example  circle  radii,  relative , RMS  error  and  mean  current 
removed  for  each  inertial  period  and  for  an  exact  circle 


INERTIAL 

“HP" 

CIRCLE 

CURRENT 

STATION 

PERIOD 

ERROR 

RADIOS 

REMOVED 

(*3> 

(111) 

(cm 

s;i) 

u 

y 

1 

1 

.3505 

.321 

-  .13 

.  04 

2 

.3775 

.631 

-  .39 

-  .56 

3 

.3810 

.839 

-5.  18 

.38 

4 

.2622 

1.  152 

1.15 

1.69 

5 

.1659 

2.663 

-6.53 

1.33 

6 

.1688 

2.602 

-  .52 

2.03 

7 

.1370 

2.859 

6.18 

-10.77 

8 

.1580 

2.577 

1.22 

2.  44 

9 

.2110 

1.128 

.24 

.  43 

10 

.2219 

1.560 

-3.23 

4.73 

1 1 

.2165 

1.714 

1.16 

-2.  62 

12 

.2695 

1.302 

2.89 

-1.64 

3 

1 

.0813 

2.547 

-  .01 

.29 

2 

.0738 

2.335 

1.83 

-  .  29 

3 

.1501 

2.  106 

.19 

-1.  17 

4 

.1530 

2.035 

.20 

-  .92 

5 

.2813 

1.440 

-5.67 

2.92 

6 

.0339 

12.455 

-9.66 

-2.  81 

7 

-0322 

12. 116 

-1.83 

5.  89 

8 

.0160 

13.394 

3.52 

4.75 

9 

.0316 

1 0.904 

3.56 

-  .16 

10 

.0301 

6.357 

2.07 

3.77 

11 

.0583 

5.818 

-  .67 

.  89 

12 

.0859 

3.891 

-2.00 

-2.  16 

True  Ci 

rcle 

.0033 

.300 

.0022 

.500 

.0002 

1.000 

.0003 

2.000 

.0002 

3.000 

.0005 

4.000 

.0023 

8.000 

.0002 

16.000 

.0021 

3  0.000 

amplitude  of  the  inertial  current  of  about  90  ca  s~l.  The 
energy  values  at  station  one  are  auch  smaller  than  those  at 
station  two  and  three.  The  technique  employed  to  fit  the 
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true  circle  to  the  manipulated  data  accounts  for  some  of 
this  difference.  When  the  data  points  are  not  circular,  the 
circle  is  adjusted  to  some  value  less  than  the  largest 
dimension  of  the  plotted  data.  When  the  data  form  a  circle, 
the  largest  dimension  of  the  plotted  data  and  the  dimension 
of  the  circle  nearly  coincide.  Por  stations  two  and  three, 
a  sharp  increase  in  kinetic  energy  is  indicated  in  Pig.  17 
during  the  fifth  inertial  period.  These  energy  values 
remain  high  for  two  to  three  iner*ial  periods  before  rapidly 
decreasing  toward  the  pre-storm  levels.  Similar  plots  for 


the 

neat 

deepest  level 

at  each  station 

are  shown  i 

Fig. 

18a. 

At  stations  one 

and  two  the  energy 

peaks  rapidly 

falls  off  rapidly  and  then  increases  to  the  end  of  the 
record.  At  station  three,  a  slower  increase  is  noticed,  and 
is  followed  by  a  nearly  constant  value  to  the  end  of  the 
record.  A  slower  increase  in  energy  is  found  (Pig.  18b)  at 
the  437  m  and  457  a  levels  at  station  three.  The  energy 
values  do  not  appear  to  have  peaked  by  the  end  of  the 
record.  Since  the  records  of  Pig.  18  are  within  20  a  of 
each  other,  it  would  be  reasonable  to  assume  the  records 
would  be  nearly  identical.  The  difference  could  be  due  to 
bottom  reflection  of  energy  since  the  457  m  instrument  was 
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positioned  only  8  aeters  froa  the  bottoa.  At  the  very 
least,  the  procedures  used  to  deteraine  the  energy  associ¬ 
ated  with  the  inertial  aotion  should  produce  results  accu¬ 
rate  to  within  a  factor  of  two  and  should  also  indicate  the 
trends.  No  confidence  is  placed  in  the  details  of  the  fluc¬ 
tuating  values.  However,  the  large  peaks  and  long  term 
increases  shown  in  the  energy  records  are  thought  to  be 
significant. 
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Pigure  17.  Energy  (a*  s-2)  versus  inerpial  period  at  the 
near  surface  levels  at  stations  one,  two  and 
three. 
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Figure  1  9. 


lanipulated  real  data  record  and  calculated 
rircle  for  (an  inertial  period  5  at  station 
ind  (b)  inertial  Dsriod  6  at  station  1. 
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A.  BACKGROUND 


Adamec  et  §1  (1981)  tested  an  embedded  aired  layer- 
ocean  circulation  model.  This  aodel  embedded  Garwood's 
(1977)  aodel  for  predicting  aired  layer  depth  and  jumps  in 
temperature  and  velocity  at  the  base  of  the  mired  layer  in  a 


six- level,  primitive  equation  o pen-ocean  model  developed  by 
Haney  (1981).  These  authors  remarked  in  their  conclusion 
that  additional  tests,  "including  comparisons  with  observa¬ 
tional  data",  were  being  planned.  This  thesis  constitutes 
an  important  step  in  those  plans  because  of  the  availability 
of  suitable  observational  data.  The  modal  used  in  the  1981 
study  is  expanded  to  simulate  the  three-dimensional  ocean 
respcnse  to  a  translating  hurricane-type  forcing  field. 
This  model  forcing  represents  an  axisy mmetric  hurricane 
translating  at  a  constant  speed  of  7.5  a  s~l.  The  model 
ocean  basin  is  a  960  km  square  with  sir  layers  of  varying 
depth  between  the  surface  and  a  free-slip  bottom  at  400  m. 
The  horizontal  resolution  is  15  km  on  a  65  by  65  grid.  The 
level  depths  and  thicknesses  are  represented  in  Pig.  21. 
The  time  step  used  in  the  1981  study,  450  s,  is  increased  to 
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600  s.  A  brief  explanation  of  the  model  as  used  in  this 
thesis  is  given  below. 

B.  MODEL  PORMOLATIOHS  A  HD  BOUNDARY  CONDITIONS 

The  governing  equations  are  the  primitive  Navier-Stokes 
equations,  the  continuity  equation,  the  hydrostatic  equation 
and  the  equation  of  state  (see  Adasec  at  al,  1981)  .  The 
hydrostatic  assumption  is  applied,  and  the  ocean  is  assumed 
incompressible  with  density  being  a  linear  function  of  temp¬ 
erature  only.  The  coriolis  force  varies  appropriate  to  the 
latitude  of  the  domain  from  25aN  to  33.6°N.  The  change  in 
coriolis  parameter  (f)  with  latitude  is  computed  as  a  finite 
difference  at  each  north-south  grid  point. 

There  are  no  fluxes  of  mass,  momentum  or  heat  across  the 
bottom  (flat)  or  side  (vertical)  boundaries.  Also,  the 
rigid  lid  approximation  (w  =  0  at  z  =  0)  is  made,  which 
requires  the  vertically  averaged  motion  be  zero.  Applying 
this  approximation  to  the  Na vie r-Stokes  equations  provides 
prediction  equations  for  the  vertical  shear  currents  (see 
Adamec  et  ^1,  1981).  Vertically  averaging  the  continuity 
equation  over  the  mixed  layer  and  applying  w  =  3  at  z  -  0 
yields  a  prognostic  equation  for  the  vertical  velocity  at 
the  base  of  the  well-mixed  layer,  w(-h).  The  entrainment 
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.  Representation  of  the  six  Bixei-layer--ocean 
circulation  uodel  layers.  Coaputation  revels 
for  u,v  and  2  ars  aiapoints  noted  in  each  level 


Figure  21 


velocity  is  predicted  by  the  aired  layer  aodel.  Together 
with  w(-h)  t  this  yields  a  prognostic  equation  for  aired 
layer  depth. 

The  axisyaaetrical  portions  of  the  hurricane  forcing  are 
shown  in  Pig.  22.  The  tangential  and  radial  stresses  are  of 
the  saae  fora.  They  are  calculated  froa  a  prescribed  wind, 
which  increases  linearly  froa  the  stora  aye  wall  to  the 
radius  of  aariaua  winds,  and  then  decreases  as  r-1/2  to  r  = 
360  ka  and  then  linearly  to  r  =  450  km.  The  radius  of  maxi- 
aua  winds  is  45  ka  and  the  inner  boundary  of  the  eye  wall  of 
the  stora  is  4.5  km.  This  wind  profile  results  in  a  wind 
stress  curl  which  is  zero  inside  the  eye  wall,  increasing 
linearly  froa  zero  at  the  eye  wall  to  a  maximum  at  45  km, 
zero  from  45  to  360  ka,  and  negative  from  360  ka  to  the 
boundary  of  the  wind  stress  (about  450  km).  The  maximum 
tangential  stress  corresponds  to  a  maximum  wind  speed  of  50 
m  s_1.  The  small  value  of  radial  stress  (-12.9  iPa  max)  is 
due  to  cross  isobaric  flow  of  about  20°. 

Eelow  the  mixed  layer,  vertical  diffusion  is  applied  to 
the  momentum  or  temperature  equations  with  a  vertical  eddy 
viscosity  or  eddy  conductivity  coefficent  (both  equal  to  0.5 
cm  s“*  at  all  depths  below  the  surface  mixed  layer).  In  the 


56 


r 


Radius  [km] 


Figure  22. 


Forcing  functions  for  the  ocean  model,  each 
normalized  by  the  value  of  the  function  at  the 
radius  of  maximum  winds.  The  tangential  and 
radial  stress  components  have  the  same  form 
(solid  line)  with  maximum  values  of  35.9  and 
-12.9  dPa.  respectively.  The  surface  heat  flux 
(dashed)  has  a  maximum  value  of  -840  W  m 
(Adamec  £t  ^i#  1981)  . 


mixed  layer,  however,  the  above  formulation  is  not 
appropriate.  The  depth  of  the  mixed  layer  is  the  boundary 


between  intense  turbulence  and  the  much  less  turbulent 


waters  beneath. 


C.  ENTRAINMENT  AND  MIXED  LAYER  MODEL  FD3MDL ATION 

In  the  mixed  layer,  prognostic  equations  for  the  mixed- 
layer  average  (bulk)  values  of  the  vertical  component  of 
turbulent  kinetic  energy  and  the  total  turbulent  kinetic 
energy  are  derived  using  the  bulk,  second-order  closure 
methods  of  Garwood  (1977).  Computation  of  the  entrainment 
buoyancy  flux  allows  calculation  of  the  downward  fluxes  of 
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heat  and  momentum  associated  with  entrainment  at  the  base  of 
the  mixed  layer. 

The  above  mixed- layer  formulation  allows  the  fraction  of 
turbulent  kinetic  energy  available  for  mixing  at  the  base  of 
the  mixed  layer  to  be  other  than  a  constant.  This  turbulent 
kinetic  energy  is  diagnostically  dependent  upon  the  surface 
buoyancy  flux  and  the  surface  friction  velocity  (known 
boundary  conditions).  It  is  thus  possible  to  couple  the 
diagnostic  mixed-layer  formulation  with  the  prognostic  ocean 
circulation  model  in  a  numerically  feasible  fashion. 

Eoth  increases  and  decreases  in  mixed  layer  depth  must 
be  considered  in  any  mixed  layer  model.  The  easier  of  these 
two  events  to  formulate  is  the  increasing  mixed  layer  depth 
case.  The  method  requires  determination  of  the  entrainment 
heat  flux  and  then  imposing  this  heat  flux  on  the  given 
temperature  profile.  Added  vertical  resolution  near  the 
base  of  the  mixed  layer  is  given  since  the  base  is  not 
required  to  coincide  with  any  of  the  prescribed  model  lev¬ 
els.  This  increased  vertical  resolution  is  very  important 
since  the  thermocline  profile  determines  the  potential 
energy  of  the  upper  ocean.  The  dynamics  of  the  mixed  layer 
and  the  ocean  circulation  are  dependent  on  this  potential 
energy. 
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Formation  of  a  new,  shallower  mixed  layer  occurs 
whenever  the  turbulent  mixing  is  not  able  to  penetrate  all 
the  way  to  the  previously  established  depth  of  the  mixed 
layer.  This  layer  reformation  occurs  when  warming  of  the 
surface  layer  occurs  with  no  increase  in  the  wind  shear  to 
sustain  the  mixing.  Numerical  formulation  of  this  event  is 
difficult  because  the  previous  structure  at  the  interface 
can  not  readily  be  preserved.  Thus  energy  budget  problems 
may  arise  when  deepening  to  the  prior  mixed  layer  depth  does 
occur.  This  model  uses  a  numerical  procedure  (see  Adamec  et 
al,  1981)  that  preserves  potential  energy  to  ensure  the 
deepening  rate  is  as  correct  as  possible  when  the  layer 
again  deepens  to  the  earlier  mixed  layer  depth.  This  fea¬ 
ture  of  the  model  is  not  tested  here  because  the  solar  flux 
is  zero  throughout  the  simulation. 

D.  DYNAMIC  STABILITY  CONDITION 

It  is  assumed  that  the  mixed  layer  is  dynamically  unsta¬ 
ble  and  the  underlying  water  column  is  normally  dynamically 
stable.  However,  dynamic  instability  of  the  underlying 
water  column  can  sometimes  occur.  In  this  model,  vertical 
fluxes  of  heat  and  momentum  between  levels  are  imposed  so 
that  the  gradient  Richardson  number  remains  greater  than  or 


equal  to  a  critical  value.  This  generalization  of  the  aore 
comacn  convective  adjustment  is  referred  to  by  Adanec  et  al 
(1981)  as  "dynamical  adjustment". 

E.  COUPLING  OF  THE  DYNAMICAL  AND  MIXING  PROCESSES  IN  THE 
MODEL 

t 

The  mixed  layer  and  the  dynamic  portion  of  the  model, 
which  is  a  level  model  that  predicts  the  average  of  a  quan¬ 
tity  in  a  layer,  are  coupled  in  two  phases.  First,  advec- 
tive  and  diffusive  changes  in  the  upper  ocean  are  calculated 
in  the  dynamic  part  and  put  into  a  form  useable  by  the  mixed 
layer  model.  Then,  the  surface  flux  and  entrainment  changes 
are  calculated  by  the  mixed  layer  model  and  transmitted  to 
the  dynamic  part  of  the  model.  A  special  treatment  of  the 
level  which  contains  the  base  of  the  mixed  layer  is  required 
(see  Adamec  et  al,  1981)  . 

F.  COMPARISON  OF  THE  MODEL  AND  THE  DATA 

The  model  provides  a  four-dimensional  space-time  view  of 
the  ocean  response  to  a  moving  hurricane.  The  data  provided 
by  NAVOCEANO  is  two-dimensional,  time  and  depth,  at  three 
different  stations.  All  three  locations  are  to  the  right  of 
the  storm  track  beyond  the  radius  of  maximum  winds  of  the 
storm.  The  differences  in  instrument  depths  and  the  limited 


number  of  stations  do  not  allow  the  extension  of  the  data 


into  a  third  dimension.  The  comparison  of  the  data  and  the 
model,  then,  is  an  attempt  to: 


•  compare  the  general  physical  characteristics  of  the 
model  results  to  see  if  they  are  reasonable  in  light  of 
the  observed  data  results. 


•  choose  a  data  set  from  the  model  which  is  similar  to 
the  raw  data  in  number,  depth  and  location  relative  to 
the  store. 


•  compare  this  data  to  the  NAVOCEANO  data  to  show  whether 
or  not  the  model  j.s  providing  a  realistic  response  to 
the  simulated  hurricane  forcing,  and 


•  examine  the  four-dimensional  model  results  for  charac¬ 
teristic  ocean  response  to  hurricane  forcing  which 
could  not  be  observed  at  the  data  stations. 


G.  GENERAL  CHARACTERISTICS  OF  THE  MODEL 


The  initial  position  of  the  hurricane  is  at  grid  point 


(23.16) .  Since  the  grid  spacing  is  15  km,  grid  point 

(23.16)  corresponds  to  a  distance  in  the  x  direction  of  33C 


km  and  in  the  y  direction  of  225  km.  Model  output  includes 


north-south  and  east-west  velocity  components,  y  and  u  ; 
temperature,  2  ;  and  mixed  layer  depth,  g.  The  initial  val¬ 
ues  of  these  quanties  are  given  in  the  following  table. 

The  quantities  a,  v,  and  J  are  extracted  at  both  three  hour 
intervals  for  the  entire  grid  at  all  six  levels  and  at  ten 
minute  intervals  for  selected  grid  points  and  levels.  The 
mixed  layer  depth  (MLD)  is  stored  at  three-hour  intervals. 


61 


TABLE  III 

Initial  values  of  modal  output  variables 

T  at  level 

HI  31  l  l  4  1  6 

(cm  s-»)  (a)  (°C) 

0  0  30.  30.0  29.2  28.  1  23.4  18.4  13.1 

The  results  predicted  by  the  nodel  at  35  and  48  h  at  level  1 
are  presented  in  Figs.  23a--d  and  Figs.  24a — d,  respec¬ 
tively.  Osing  the  storm  translation  speed  of  7.5  a  s-1,  the 
distance  of  the  storm  from  the  bottom  of  the  grid  is  given 
in  Table  IT.  The  storm  center  has  tracked  off  the  grid  by 
hour  24  and  the  southern  boundary  of  the  storm  has  passed 
the  edge  by  hour  48. 

TABLE  IT 

Hurricane  distance  from  bottom  of  grid  at  six-hour  intervals 


DISTANCE 

HOUR 

330 

0 

492 

6 

654 

12 

816 

18 

978 

24 

1140 

30 

1302 

36 

1464 

42 

Several  important  features  are  noticeable  in  Figs.  23 


and  24: 


•  The  current  fields  at  36  and  48  hours  are  "out  of 
phase”,  i-3-  the  areas  of,  westerly  currents  at  hour  36 
are  generally  easterly  at  hour  4  8. 
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TSPO  7.5  M/S  TSPO  7.S  M/S 

U  LEVEL  1  H8UR  36  V  H0uR  36 


Figure  23.  Model  output  of  (a)  u  velpcity  coaponent,  (b)  v 
velocity  coaponent,  (c)  aixad-layer  depth,  and 
(d)  layer  teaperature  at  36  hours. 
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Figure  24. 


Hodel  output  of  (a)  u  velocity  conponent,  (b)  \ 
velocity  component,  (c)  aixed-layer  deoth,  and 
(d)  layer  temperature  at  48  hours. 
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•  Current  velocities  in  the  aixed  layer  range  froa  neg¬ 
ative  values  of  about  70  ca  -l  to  positive  values  of 
about  80  ca  s~*. 


•  The  aired  layer  deepens  preferentially  to  the  right  of 
the  storm,  and  ranges  froa  ainiaum  values  of  about  30  a 
to  aaximua  values  of  alaost  60  a. 


•  A  ridge  of  cold,  upwelled  water  lies  m  the  wake  of  the 
store.  This  ridge  is  shown  best  in  the  level  3  temper¬ 
ature  contour  of  Pig.  25.  Superposed  on  this  cold 
ridge  is  an  upwelling /downwelling  wake  with  an  inertial 
period,  Opwelling  is  in  the  areas  of  ainiaum  MLD's 
(32.5  m  Pig.  24c),  and  downwelling  is  in  the  area  of 
aaxiaua  MLD's  (58.0  in  Pig.  23c).  Minimum  MLD's  are 
reached  at  the  end  of  the  upwelliag  cycle  and  maximum 
MLD's  are  reached  at  tha  end  of  the  downwelling  cycle. 


These  results  are  consistent  with  the  results  derived  from 
the  NAVOCSANO  data,  which  shows  that  the  primary  ocean 
current  response  to  hurricane  passage  whan  no  mean  current 


exists  is  circular  motion  at  tha  local  inertial  period.  The 
reversal  in  flow  direction  every  12  hours  is  particularly 
evident  in  the  u  velocity  component  records  for  level  six 
shown  in  Pig.  26.  The  model  predictions  describe  the  ocean 


response  over  the  entire  grid,  rather  than  just  at  some 
selected  points.  Trajectories  are  plotted  at  several  loca¬ 
tions  perpendicular  to  the  storm  track  (see  Fig.  27)  .  The 


surf  ace- layer  trajectory  that  began  on  the  left  side  of  the 
track  (Pigs.  27a)  clearly  shows  the  surge  to  the  left  as 
the  storm  approaches,  strong  inertial  aotion  following  storm 


passage  and  a  net  transport  to  the  south.  The  trajectory 
along  the  storn  path  (Pig.  27b)  shows  a  surge  to  the  left, 
recovery  toward  the  right  and  a  snail  net  displacement.  The 
trajectories  that  began  on  the  right  (Pigs.  27c  and  27d) 
show  a  net  deflection  northward  and  to  the  right.  These 
trajectories  show  that  different  transport  profiles  are 
expected  at  different  locations  relative  to  the  storm.  This 
feature  is  not  obvious  in  the  N AVOCE AND  data  due  to  the 
limited  number  of  locations  of  the  observations. 

H.  SELSCTED  POSITION  AND  DEPTH  DATA 

A  data  set  similar  to  the  NA  70CEAN0  data  set  was  derived 
by  extracting  model  variables  at  ten  minute  intervals  at  the 
three  positions  and  depths  shown  in  Table  V.  These  time 
series  of  model  variables  (Pigs.  28  and  29)  are  plotted  sim¬ 
ilarly  to  the  raw  data  plots  of  Figs.  6  and  7.  The  observed 
(Pigs.  6)  and  predicted  (Fig.  28)  u-ooiponents  show  very 
similar  variations  although  the  model  data  are  much 
smoother.  Also,  although  the  magnitude  of  the  velocity  com¬ 
ponent  at  level  2  of  Pig.  28a  is  somewhat  smaller  than  the 
21  a  magnitude  shown  in  Pig.  5a,  the  magnitude  at  level  6 
(300  a)  of  Pig.  28b  is  much  smaller  than  the  U37  m  magni¬ 
tude  shown  in  Pig.  6b.  The  implication  is  that  not  enough 
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Figure  27 


energy  is  being  vertically  transmitted  by  the  aodel.  Both 
the  aodel  aotion  and  the  observed  floe  are  inertial  and  nei¬ 


ther  has  completely  daaped  by  the  end  of  nine  days.  The 
model  u  components  of  the  velocity  at  levels  one  and  six  are 
exactly  180  degrees  out  of  phase  (Pig.  30a)  .  That  is^,  the 
velocities  of  the  surface  and  bottom  model  layers  are  in  the 
opposite  directions.  The  NAVOCEANO  velocities  are  also 
exactly  out  of  phase  immediately  after  storm  passage, 
although  they  become  more  in  phase  later  in  the  record  (see 
Pig.  30b).  For  any  given  level  of  the  aodel,  the  u  and  v 
variations  are  90  degrees  out  of  phase. 

TABLE  V 

Positions  and  depths  for  model  ten  minute  interval  data  and 
corresponding  NAVOCEANO  data. 

DISTANCE 


INERTIAL 

FROM 

MID 

NAVOCEANO 

GRID 

PERIOD 

STORM  CENTER 

LAYER 

INSTRUMENT 

POSITION 

MODEL  NAVOCEANO 

MODEL  NAVOCEANO 

DEPTH 

DEPTH 

(*#y) 

(hours) 

(km) 

(m) 

(a) 

30,33 

24.50  24.54 

105  100 

37.5 

49. 

75.0 

64. 

150.0 

92. 

2  9,35 

24.30  24.34 

90  85 

19.0 

19. 

150.0 

179. 

300.0 

324. 

31,39 

23.95  23.90 

120  125 

19.0 

21. 

75.  0 

251. 

150.0 

437. 

300.  0 

457. 

Three 

-hour  running  averages  by  inertial  period 

were  com- 

puted  for 

the  extracted 

aodel  data  using 

the  same 

technique 
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Figure  28. 


0  coipocents  of  the  aodel  current  output  for  1C 
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Figure  30.  Superpositions  of  the  surface  and  deep  layer 
plots  xn  (a)  Fig.  28  and  (b)  Fig.  6. 


as  used  for  the  NAVOCEANO  results.  A  direct  comparison  of 
Figs.  13,  14  and  15  to  the  extracted  model  results  running 
inertial  period  averages  is  made  in  Figs.  31,  32,  33.  The 
nodel  near-surface  residual  flows  are  very  similar  to  each 
other  (Figs.  31a  and  32a)  ,  as  are  the  lower  level  results 
(Figs.  31b,  32b  and  33).  Net  transport  to  the  north  is 
indicated  in  all  records  of  the  model.  The  NAVOCEANO  data 
indicate  net  transports  that  are  closely  tied  to  the  local 
topography,  which  is  absent  in  the  model.  The  NAVOCEANO 
results  also  show  opposite  current  directions  for  near-sur¬ 
face  and  aid  depth.  These  observations  provide  interesting 
examples  of  results  not  shown  in  the  model,  and  which 
require  addtional  investigation  to  explain. 

In  order  to  make  a  more  quantitative  comparison  between 
the  model  and  the  NAVOCEANO  inertial  oscillations,  the  model 
results  from  the  selected  stations  shown  in  Table  V  were 
treated  in  the  same  manner  as  the  NAVOCEANO  raw  data.  Cir¬ 
cle  radii,  relative  RMS  error  and  mean  current  removed  are 
listed  (Table  VII)  similarly  to  Table  II.  Relative  RMS 
errors  are  similar  to  those  for  station  three  (Table  II) 
although  maximum  circle  radii  are  not  as  large.  Resultant 
calculations  of  energy  (cm2  s-2)  versus  inertial  period  are 
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plotted  (Figs.  34  and  35)  similarly  to  the  energy  plots  of 
Figs.  17  and  18.  These  figures  show  that  the  major  differ¬ 
ences  in  the  aodel  and  the  NAVOCEANO  data,  at  least  in  a 
qualitative  sense  are  in  the  rate  of  propagation  of  energy 
with  depth.  It  appears  that  less  energy  is  propagated  ver¬ 
tically  in  the  model  compared  to  the  observations.  The 
energies  in  Figs.  35a  and  35b  are  much  smaller  in  comparison 
to  the  values  in  Fig.  34  than  are  the  similar  results  of  the 
observed  data.  Also,  there  is  no  slow  increase  in  the  hori¬ 
zontal  kinetic  energy  near  the  bottom  in  the  model  as 
observed  in  the  two  deepest  NAVOCEANO  records.  The  level 
four  records  in  Fig.  34a  suggest  that  energy  is  being  con¬ 
tinually  propagated  into  the  layer  because  the  energy  levels 
do  net  appreciably  decay.  Detailed  calculations  are  neces¬ 
sary  to  establish  the  energy  flux  mechanisms  from  the  model 
data. 
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TABLE  VI 


Saaple  of  circle  radii,  relative  BHS  error  and  lean  current 
reaoved  for  each  inertial  period  from  the  nodel  siaulation. 


lBP 

ST  AT  10  H 
(1, 1) 


INERTIAL  “W ,E  CIRCLE 
PERIOD  ERROR  RADIOS 

<*■)  <*a> 


C05FSNT 
BEHOVED 
(£t  Sli) 


(30,33) 


1 

.0580 

6.373 

2.33 

2.73 

2 

.0334 

10.166 

-5.28 

5.63 

3 

.0430 

6.877 

-1.86 

-  .  14 

4 

.0532 

5.474 

-1.53 

.08 

5 

.0495 

4.817 

-  .93 

.  25 

6 

.0621 

4.  173 

-  .76 

.33 

7 

.0604 

3.650 

-  .67 

.  35 

8 

.0510 

3.187 

-  .63 

.  31 

9 

.0679 
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-  .36 
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Figure  34. 


Energy  (a*  s-*)  versus 
near  surface  levels  oi 
and  31,39. 
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IV.  SOHMARX.  CONCLUSIONS  MJ>  RECOMMENDATIONS 


A.  SUMMARY 

The  NAVOCEANO  data  provided  several  physically 
significant  results.  The  currents  measured  by  the 
instruments  were  affected  by  topography  and  by  the  storm. 
The  flow  tended  to  be  along  the  bathymetry  contours  except 
that  the  storm  imparted  a  surge  to  the  flow  consistent  with 
wind  stress  associated  with  the  storm.  The  station  one 
records  do  not  depict  easily  identifiable  inertial  motion. 
Station  one  is  in  shallow  water  of  about  100  meters,  and 
shoreward  motion  to  mid-depth  is  evident  prior  to  the  storm 
passage.  An  offshore  flow  is  evident  from  about  the  time  of 
the  storm  passage  for  a  duration  of  about  two  inertial  peri¬ 
ods.  The  tendency  at  each  station  is  for  the  post-storm 
flow  to  return  to  the  same  direction  as  the  pre-storm  flow. 
Energy  input  by  the  storm  has  generally  dissipated  by  a  fac¬ 
tor  of  e-»  after  three  or  four  inertial  periods,  and  the 
post-storm  flow  magnitude  is  greater  than  the  pre-storm  mag¬ 
nitude.  The  horizontal  kinetic  energy  associated  with  iner¬ 
tial  motion  in  the  near-surface  water  increased  more  than 
two  orders  of  magnitude  as  the  storm  passed  and  then  decayed 
over  six  or  seven  inertial  periods.  At  the  deepest 
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instrument  levels,  the  kinetic  energy  associated  with  iner¬ 
tial  motion  increased  slowly  and  had  not  peaked  seven  iner¬ 
tial  periods  after  storm  passage.  The  mid-level  instruments 
recorded  either  a  sharp  peak  after  storm  passage  or  a  less 
rapid  increase  and  little  decay  to  the  end  of  the  record 
depending  upon  location.  The  near-surface  and  near-bottom 
flows  are  about  180°  out  of  phase  for  at  least  five  or  six 
inertial  periods  after  storm  passage. 

The  mixed  layer-ocean  circulation  model  results  are  very 
similar  in  many  respects  to  the  NA7OCEAN0  observations.  The 
model  does  not  include  bathymetry,  and  therefore  bathymetric 
effects  are  not  reproduced.  The  model  does  show,  however,  a 
storm-induced  surge,  inertial  motion  damping  over  several 
inertial  periods  at  the  surface  and  bottom  layers  and  dif¬ 
ferent  magnitudes  in  pre- storm  and  post-storm  flow.  The 
model  predictions  indicate  an  almost  instantaneous  near-in¬ 
ertial  energy  propagation  rats  with  depth.  The  flow  at  the 
surface  and  at  depth  have  the  same  form,  increasing  in  mag¬ 
nitude  and  decaying  at  very  nearly  the  same  rate  at  each 
level.  At  mid  levels  of  the  model,  the  kinetic  energy  level 
remains  nearly  constant.  The  current  oscillations  in  the 
mixed  layer  and  in  the  bottom  layer  are  180°  out  of  phase. 


Trajectories  of  the  model  resalts  show  that  the  net  dis¬ 
placement  at  some  stations  is  nearly  zero  while  at  other 
stations  the  net  displacement  is  tens  of  kilometers.  Also, 
the  direction  of  surge  and  of  the  residual  flow  is  dependent 
on  location  relative  to  the  star  a  track.  The  model  shows  a 
zone  of  upwelling/dovnwelling  in  the  wake  of  the  storm  which 
has  a  near-inertial  frequency. 

B.  CONCLUSIONS 

The  bathymetry  of  the  ocean  is  very  important  in  determ- 
ing  the  direction  of  flow  of  currents  at  all  levels.  This 
part  of  the  flow  is  not  included  in  this  ocean  circulation 
model  and  must  be  subtracted  from  the  real  ocean  data  if 
consistent  comparisons  of  model  results  and  current  observa¬ 
tions  are  to  be  made.  The  effect  of  the  storm  is  to  drive 
the  average  motion  in  a  direction  consistent  with  the  wind 
field  of  the  storm.  Changes  in  direction  occur  over  a  very 
short  time  period  and  are  observed  at  all  depths.  As  the 
storm  passes,  inertial  motion  results  due  to  a  balance 
between  the  centrifugal  and  coriolis  accelerations.  The 
kinetic  energy  associated  with  this  motion  remains  nearly 
constant  near  the  surface  for  two  to  three  inertial  periods 
and  decays  over  an  additional  three  to  four  inertial  periods 
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to  a  post-store  level  greater  than  the  pre-store  level.  The 
available  data  do  not  show  the  energy  levels  returning  to 


the  Fre- store  level  before  the  end  of  the  record.  There  is 
an  obvious  tiee  lag  of  the  store  which  increases  with  depth. 
The  exact  tiee  of  the  initial  effect  at  any  particular  depth 
is  difficult  to  measure  since  the  magnitude  of  the  velocity 
does  not  instantaneously  increase,  but  increases  slowly  with 
the  onset  of  the  storm.  Nonetheless,  tha  downward  propaga¬ 
tion  of  this  energy  is  net  instantaneous  but  requires  on  the 
order  of  one-half  inertial  period  to  propagate  one-half  kil¬ 
ometer.  The  energy  propagates  to  a  depth  of  about  450  m  in 
one-half  of  an  inertial  period.  At  this  depth  tha  energy 
increases  for  at  least  seven  to  eight  inertial  periods.  The 
mixed  layer-ocean  circulation  model  closely  simulates  the 
inertial  motion  and  damping  rates  near  the  surface.  The 
vertical  propagation  of  energy  is  not,  however,  very  similar 
to  the  NAVOCEANO  observations.  The  model  does  provide 
insight  into  th“  three-dimensional  ocean  and  is  impressive 
in  distinguishing  between  results  at  locations  varying  with 
respect  to  the  store  center. 
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C.  BSCO  MMEN  DATION S 


The  primary  differences  between  the  MAVOCEANQ  data  and 
the  aodel  results  are  in  the  vertical  propagation  of  kinetic 
energy  associated  with  storm  passage.  Further  research  is 
required  to  isolate  the  mechanism  of  energy  propagation  so 
that  it  can  be  included  in  the  aodel  formulation. 
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